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Abstract 

The effects of both eUiptical shape and stage of emergence of the coronal loop on the 
resonant absorption of standing kink oscillations are studied. To do so, a typical coronal 
loop is modeled as a zero-beta longitudinally stratified cylindrical magnetic flux tube. We 
developed the connection formulae for the resonant absorption of standing transversal oscil- 
lations of a coronal loop with an elliptical shape, at various stages of its emergence. Using 
the connection formulae, the dispersion relation is derived and solved numerically to obtain 
the frequencies and damping rates of the fundamental and first-overtone kink modes. Our 
numerical results show that both the elliptical shape and stage of emergence of the loop 
alter the frequencies and damping rates of the tube as well as the ratio of frequencies of the 
fundamental and its first-overtone modes. However, the ratio of the oscillation frequency to 
the damping rate is not affected by the tube shape and stage of its emergence and also is 
independent of the density stratification parameter. 
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1 Introduction 



Solar corona and its extraordinary high temperature has been the topic of various debates and 
studies from several decades ago. The origin and the source of coronal continual heating and 
high temperature have been related to coronal loops. The claim that the coronal loops and their 
behaviors such as their damping oscillations may be one of the main reasons of coronal heating, 
has been investigated through several studies so far. 

Transverse oscillations of coronal loops have been observed by the Transition Region and 
Coronal Explorer (TRACE) for several years (see e.g. Aschwanden et al. 1999; Schrijver &; 
Brown 2000). Nakariakov ct al. (1999) interpreted these oscillations as fast kink modes with the 
period ranging from 2.3 to 10.8 min and decay time from 3.2 to 20.8 min. The observed values 
of the periods and decay times make it possible to obtain indirect information on the conditions 
of the plasma and magnetic field in coronal loops. 

Ofman &: Aschwanden (2002) used the data deduced by Aschwanden et al. (2002) to inves- 
tigate the oscillations of 11 coronal loops. They argued that the observed TRACE loops consist 
of multiple unresolved thin loop threads which produce inhomogeneous internal structure of the 
observed loop. They adopted 1-dimensional cartesian slabs of plasma with the magnetic field 
lines in the z-direction and the direction of the inhomogeneity along the x-axis normal to the 
magnetic surfaces, as a simple model for the oscillating loops. They found that the dependence 
of the decay time on both the length L and the width w of the loop is in excellent agreement 
with the power law damping predicted by phase mixing. 

The property of resonant absorption as a non-thermal mechanism makes it possible to de- 
scribe the heating of magnetic loops in solar corona as well as rapid decaying of magnetohydro- 
dynamics (MHD) waves even in weakly dissipative plasmas (see e.g. lonson 1978; Poedts et al. 
1989; Ofman et al. 1994; Erdclyi & Goossens 1994, 1995; Tirry & Goossens 1996; Andries et al. 
2005b; Safari et al. 2006; Dymova & Ruderman 2006; Goossens et al. 2009). 

Verwichte et al. (2004), using the observations of TRACE, have identified the fundamental 
and its first harmonic of the transverse kink mode in two coronal loops. The period ratios 
observed by Verwichte et al. (2004) are 1.81 =b 0.25 and 1.64 it 0.23. However, these values were 
corrected with the improvement of the observational error bars to 1.82 it 0.08 and 1.58 ± 0.06, 
respectively, by Van Doorsselaere et al. (2007). Also Verth et al. (2008) added some further 
corrections by considering the effects of loop expansion and estimated a period ratio of 1.54. All 
these values clearly are lower than 2. This may be caused by different factors such as the effects 
of density stratification (see e.g. Andries et al. 2005a; Erdelyi k. Verth 2007; Karami & Asvar 
2007; Safari et al. 2007; Karami et al. 2009) and magnetic twist (see Erdelyi & Carter 2006; 
Erdelyi & Fedun 2006; Karami k Barin 2009; Karami & Bahari 2010, 2012) in the loops. 

Karami et al. (2009, hereafter Paper I) investigated the effect of longitudinally stratification 
on resonant absorption of MHD waves for both kink (m = 1) and fluting (m = 2) modes. They 
found that the frequencies and damping rates of both the fundamental and first-overtone modes 
increase when the stratification parameter increases. Also for stratified loops they obtained the 
ratio of the frequencies 0J2/ijJi of the first overtone and its fundamental mode less than 2. 

Morton k, Erdelyi (2009) studied the effects of both the elliptical shape and stage of emer- 
gence of the loops on the period ratio P1/P2 for the minor and major elliptical cases. Their 
results showed that the parameter characterising the stage of emergence does affect the value of 
period ratio Pi/P2- Particularly, the greatest contribution from emergence to the period ratio 
occurs when the loop is fully emerged. Also they showed that the ellipticity of the loop has an 
important role in the value of P1/P2 for minor elliptical case but the major ellipse was found to 
have a less effect on the period ratio of standing oscillations. 
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Here we combine the two models considered in Paper I and Morton & Erdelyi (2009) to 
investigate of the effects of both eUiptical shape and stage of emergence of the coronal loop on 
the resonant absorption of standing transversal kink oscillations observed by the TRACE. This 
paper is organized as follows. In Sections 2 and 3, we combine the two techniques of Paper I and 
Morton & Erdelyi (2009) to derive the equations of motion, introduce the relevant connection 
formulae and obtain the dispersion relation. In Section 4, we give numerical results. Section 5 
is devoted to conclusions. 



2 Equations of motion and modeling of the flux tube 

The linearized MHD equations for a zero-beta plasma are given by 

^ = 7^{(V X <5B) X B + (V X B) X 6B} + ^V2<5v, (1) 
at Anp p 



d5B 



c 



,2 



^ Vx(*vxB) + _V^ffl, (2) 

where Sv and SB are the Eulerian perturbations of velocity and magnetic fields; B, p, a, r} and 
c are the background magnetic filed, the mass density, the electrical conductivity, the viscosity 
and the speed of light, respectively. 

The simplifying assumptions are the same as in Karami &i Asvar (2007). According to 
Andries et al. (2005b) and Paper I, one can expand the perturbed quantities d\ and SB as 
follows 

oo 

6B{r,z) = ^5B('=)(r)V'(*=)(z), 
fe=i 

oo 

(5v(r,z) = ^(5v('=)(r)^/;(^)(z), (3) 
k=l 

where ip^''\z)s form a complete set of orthonormal eigenfunctions and satisfy the eigenvalue 
relation 

W=) = av'^'=\ (4) 

where La is the Alfven operator, 

with Alfven velocity va = and straight constant background magnetic filed B = Bz. 

We further assume there is a density stratification along the tube axis in ^;-direction.Since 
we are interested in resonantly damped oscillations, it implies that the density profile must be 
radially structured too. Following Paper I, we consider the density profile given by 

p{r,z) = po{r)p{z), (6) 

where 

Pin, {r < Rl), 

(i?-r) + /9ex, iRi<r<R), (7) 
Pex, {r > R). 



Po{r) 



Pin— Pex 

R-Ri 
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Here, R is the loop radius and Ri < R is the radius of the homogeneous part of the tube. 
The radius at which resonant absorption occurs is between Ri and R. The thickness of the 
inhomogeneous layer, I = R — Ri, will be assumed to be small. Here, and are the 
footpoint densities of the interior and exterior regions of the tube, respectively. 

According to Morton & Erdelyi (2009) we consider two types of elliptical loop that can occur, 
the minor elliptical loop where minor axis of the ellipse is the vertical axis of the loop, and the 
major elliptical loop where major axis of the ellipse is the vertical axis of the loop. Note that 
the minor ellipse is a situation that occurs most plausibly under coronal conditions. 

For the minor elliptical case, the longitudinally stratified density profile takes the form 



exp 



cos 



-/X- 



{a{z)) 



1 - sin^ 



A 



V 



(8) 



and for the major one, it is given by 



p{z) = exp 



cos a(z) 1 — e 



1/2 



1 - cos^ 



(a(x 



-A 



V 



where 



1-A 



,2\ V2 



a2 J 



(9) 



(10) 



is the ellipticity of the loop with minor half-axis of length b, and major half-axis of length a. 
Also n := ^ is defined as stratification parameter, where H and L are the density scale height 
and length of the loop, respectively. The parameter A describes the stage of emergence of the 
loop from the photosphere. It is defined as the ratio of the distance of the photosphere from 
center of the ellipse to the vertical half-axis. A positive value of A refers to the situation in 
which the center of ellipse is sitting below the photosphere (early stage emergence), and thus, 
the negative A for the center above the photosphere (late stage emergence). A zero value of A 
corresponds to a loop having a semi-elliptical shape. For the minor elliptical case, A is given by 



A 



(11) 



where fj,H is the distance of the loop apex from the photosphere. For the major one, A is defined 
as 



A = 1 - 



(12) 



Note we have considered a tube of length L which its footpoints are in the two points z = and 
z = L, and also note that in Eqs. (8) and (9), a{z) is the angle between the vertical axis of the 
loop and the line joining the center of the ellipse to the plasma element located at distance z 
along the tube. Following Morton Sz Erdelyi (2009) for the minor elliptical case, one can obtain 
the value of a{z) by calculating the ellipse arc length defined as 



J^' (l - £2 sin2(i)) = ( § - - 'di, 



(13) 



4 



where ti and t2 are parametric angles given by 

ti = arctan ^— tan(a)^ , t2 = arctan ^— tan(^) 



and 



9 = arctan 



1 / 1 - A 



1 ■ 

2\ 2 



A \ 1 - 



(14) 



(15) 



is the angle between the vertical axis of the loop and a line that joins the ellipse center to the 
loop foot-point (see Fig. 2 in Morton Sz Erdelyi 2009). 
For the major elliptical case, we have 



where 



and 



ti = arctan ^- cot(a)^ , t2 = arctan ^- cot(^)^ , 

((l-a(l-A^ 



6 = arctan 



1 



A 



(16) 



(17) 



(18) 



According to Paper I, in the absence of dissipation, in the interior region (r < i?i), solutions of 
Eqs. (1) and (2) are 



+ 00 



<5i?f )(r,z) = E A('°''^)Jm(|/cin,k|r-)V'^'"''^(^), 



(19) 



k=l 



6vl}-\r, z) = ?^^^^"''^^m(IMn,k|r)^(''^''^)(.), (20) 



where 



Here Jm is the Bcsscl function of the first kind and a prime on Jm and hereafter on each function 
indicates a derivative with respect to their appropriate arguments. The solutions for the exterior 
region r > R, are the same as equation (20) except that Jm, index "in", and |A;in,k| are replaced 
by Kja, "ex" and A;ex,k = — ^^4-^ ' respectively, everywhere. Where is the modified Bessel 
function of the second kind and shows that the wave amplitude vanishes in large distance away 
from the tube boundary. 
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3 Boundary conditions, Connection formulae and dispersion re- 
lation 



In the absence of dissipation effects, an appropriate dispersion relation is obtained by requiring 
that the solutions for perturbed quantities are continues at the tube surface. When a dissipative 
layer is considered, the solutions may experience jumps across the layer. An appropriate relation 
connecting the solutions of outside and inside the tube, is called the "connection formulae". 
Following Paper I, the jump across the boundary (resonance layer) for SB^ and Svr is 

[SB,] = 0, (22) 



where 



and 



with 



[6vr 



Lai = 



E 



BLbm? 




SBt'""^) 




in,k) 


Lai 


0(m,k)\> 



(in.k) 



dLA 



dr 



r=rA 



(in,k) 



, (in,k) _ 



,(in,k) 



kk\ 



sm 



\dpo{r) 



dr 



JTT 



j = k 



r=rA 



(23) 



(24) 



(25) 



(26) 



Here satisfies La4>''''^'^^ = and Cj = u + where 7 is damping rate. Also Ri < ta < R 

Ri < rA < R is the radius at which the singularity occurs. Note that I = R — Ri is the thickness 
of the inhomogeneous layer and Davila (1987) showed that in the resonance absorption, however, 
the damping rate is independent of the dissipation coefHcient values. But the resonance layer 
width scales as (5a oc (^ + ^)^^^- Karami k Bahari (2010) showed that for the Reynolds 

7^ = = 560 and Lundquist S = (^^)/(^) = 10^ numbers given by Ofman 

et al. (1994), and taking L = 10^ km, R/L = 0.01, a/R = 0.08, Pcx/Pin = 0.1, and interior 
Alfven velocity va^^ = 2000 km s~^ for a typical coronal loop, then one can get 5\ ~ 85 km 
which is very close to the thickness of the inhomogeneous layer I ~ 87 km. This suggests that 
one can use the thin boundary approximation which assumes that the thickness of the resonance 
layer is the same as the inhomogeneous layer width (see also Goossens et al. 2009). 

Note that the effects of elliptical shape and stage of emergence of the loop on the resonant 
absorption appear in jump conditions via the function Skj which is obtained as 

= "II Jo (t^) i^^'^) {t') ^^^^ 

Substituting the fields of equation (19) and (20) in jump conditions (22) and (23) gives 



/ n 



(ex,l) 



-n 



(in,l) 



n 



^(ex,l) 

(ex,l) 
2 

^(ex,l) 
"2 



"1 



(in,l) 



n 



-n 



(in,l) 



"2 



(ex,2) 
^(ex,2) 

„(ex,2) 

^(ex,2) 
"2 



_n(in.2) 
2(in,2) _|_ p(in,2) 



2(i'^'2) _^ p(in,2) 



\ 



^(in,l) 

^(cx,2) 



V 



0, 



(28) 
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where all the definitions in Eq. (28) are as those of Paper I. Note that the dispersion relation is 
obtained by requiring that the system (28) has non-trivial solutions, i.e. its determinant is zero. 

In the next section, we solve the dispersion relation numerically to obtain the frequencies oj and 
damping rates 7 of the resonantly damped MHD kink oscillations of longitudinally stratified 
elliptical emerging coronal loops. 

4 Numerical Results 

To solve the dispersion relation (28) numerically, we chose the physical parameters L = 10^ km, 
R/L = 0.01, l/R = 0.02, B = 100 G, pin = 2 x lO^^'' g cm^^ and pex/pin = 0.1. For such a loop 
one finds 1;^;^ = ^^p. = 2 x 10^ km s~^ and oja^^ '■= = 0.02 rad s~^. In what follows, we 
illustrate our numerical studies in the three separate equilibrium cases containing (i) circle-arc 
emerged loop (ii) minor elliptical semi-emerged loop and (iii) minor and major elliptical loops. 

4.1 Circle-arc emerged loop (e = 0, A 7^ 0) 

The effect of stage of emergence of the tube on both the frequencies lo and damping rates 7 are 
calculated by numerical solution of the dispersion relation, i.e. Eq. (28). In Figs. 1 and 2, the 
frequencies, damping rates and their ratio for the fundamental and first-overtone kink {m = 1) 
modes are plotted versus the stratification parameter /i for a circle-arc flux tube (e = 0) at various 
stages of emergence containing early stage emergence (A = 0.75), semi-emerged ellipse (A = 0) 
and late stage emergence (A = —0.75). Figures 1 and 2 show that (i) for a given loop shape 
parameter A, both frequencies cji, UJ2 and their corresponding damping rates |7i|, I72I increase 
when the stratification parameter /x increases, (ii) For a given fi, both frequencies and damping 
rates increase when A increases. For instance, for /j, = 0.9, the early stage emergence A = 0.75 in 
comparison with late stage emergence A = —0.75 would cause oji,uj2, |7i| and I72I to increase by 
about 10.2%, 11.7%, 10% and 14.5%, respectively, (iii) The ratio of the oscillation frequency to 
the damping rate, i^/|7|, is independent of stratification. Also the stage of emergence of the loop 
does not affect w/|7|. This means that the ratio w/|7| is independent of how the loop center is 
situated with respect to the photosphere. 

In Fig. 3, the ratio of the frequencies OJ2/001 of the first overtone and its fundamental mode is 
plotted versus the stratification parameter. Figure 3 shows that (i) for a given A, the frequency 
ratio decreases from 2 (for an unstratified loop) with increasing density stratification, (ii) For a 
given 12, the frequency ratio decreases when A decreases. The results of UJ2/001 are in agreement 
with those obtained by Morton & Erdelyi (2009). 

4.2 Minor elliptical semi-emerged loop (e 7^ 0, A = 0) 

Figures 4 and 5 as Figs. 1 and 2 display the result of frequency, damping rate and ratio oj/IjI 
but for a minor semi-emerged loop (A = 0) with different ellipticity parameters e = 0.0, 0.4 and 
0.6. Figures 4 and 5 show that (i) for a given ellipticity parameter e, both frequencies coi, 002 
and their corresponding damping rates |7i|, I72I increase when the stratification parameter n 
increases, (ii) For a given /i, both frequencies and damping rates increase with increasing e. For 
instance, for fi = 0.9 considering a loop with e = 0.6 would cause to increase wi, UJ2, |7i| and I72I 
up to 3.5%, 4%, 3.5% and 5.3%, respectively, in comparison with a semi-circular loop (e = 0). 
(iii) The ratio oo/\^\ remains unchanged by increasing the stratification parameter. Also this 
ratio is independent of the elliptical shape of the loop. 
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Figure 6 illustrates the frequency ratio 0J2/0J1 as a function of stratification parameter for a 
semi-emerged loop (A = 0) with different ellipticity parameters. Figure 6 presents that (i) for 
a given ellipticity parameter, the ratio 002/^1 decreases from 2 by increasing the stratification 
parameter and approaches below 1.5. (ii) For a given /i, the frequency ratio U2/1JJ1 decreases 
when e decreases. This is in good concord with the result of Morton &; Erdelyi (2009). 

4.3 Minor and major elliptical emerged loops (e 7^ 0, A 7^ 0) 

Figures 7, 8 and 9 make it possible to do some comparisons between major and minor elliptical 
cases. Figures 7 and 8 show that: (i) the frequencies and damping rates of both fundamental 
and first-overtone kink (m = 1) modes of a minor elliptical case are greater than those of a 
major elliptical case and also of a circular-arc semi-emerged loop (e = A = 0) studied in Paper 
I. For instance, for jjl = 0.9 for a minor elliptical loop with e = 0.6 and A = 0.75, the values of 
^1) ^2, ItiI and I72I arc 19.1%, 26.2%, 18.4% and 22.4% greater than those of a major elliptical 
case, respectively, (ii) The frequencies and damping rates of a major elliptical loop are slightly 
greater than the circular-arc semi-emerged loop, (iii) The ratio w/|7| of both fundamental and 
first-overtone modes do not affected by stratification for minor /major elliptical and circular-arc 
semi-emerged loops. 

In Fig. 9, the frequency ratio UJ2/001 versus the stratification parameter is plotted for major 
and minor cases as well as a circular-arc semi-emerged loop. Figure 9 shows that: (i) for a given 
jjL, the frequency ratio ^2/^1 of the minor elliptical loop is greater than the major one. For 
instance, for jjL = 0.9, the ratio 002/1^1 of a minor elliptical case is 5.9% greater than that of a 
major elliptical one. (ii) The ratio UJ2/1JO1 in a major elliptical loop is slightly greater than the 
circular-arc semi-emerged loop. Therefore, the major ellipse has a lesser effect on the frequency 
ratio of standing kink oscillations. This is in good accordance with that obtained by Morton &; 
Erdelyi (2009). 

5 Conclusions 

Here we investigated the effects of both ellipticity and stage of emergence on the resonant 
absorption of standing kink waves in longitudinally stratified coronal loops. We considered a 
coronal loop as a pressureless cylindrical flux tube embedded in a straight magnetic field that 
undergoes a longitudinal density stratification and radial density structuring. We extended 
the relevant connection formulae for the resonant absorption of transverse kink oscillations of a 
coronal loop with an elliptical shape, at various stages of its emergence from the sub-photosphere 
into the solar corona. We studied three stages of emergence of the loop at which the center of 
eclipse is sitting below (A > 0), on (A = 0), and above (A < 0) the photosphere. We considered 
two types of elliptical loop that can occur, the minor ellipse and the major ellipse. Note that the 
minor ellipse is a situation that occurs most plausibly under coronal conditions. By numerically 
solving the dispersion relation, we obtained the frequencies and damping rates of the fundamental 
and first-overtone kink (m = 1) modes. Our numerical results show the following. 

i) By increasing the density stratification parameter fi in the loop, both frequencies and 
damping rates increase while the frequency ratio uj2/^i decreases. 

ii) In a circle-arc emerged loop (e = 0, A 7^ 0) for a given ^u, the frequencies, damping rates 
and the frequency ratio increase when the stage of emergence parameter A increases. 

iii) In a minor elliptical semi-emerged loop (e 7^ 0, A = 0) for a given 11, the frequencies, 
damping rates and the frequency ratio increase when the ellipticity parameter e increases. 
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iv) In a minor elliptical emerged loop (e 7^ 0, A 7^ 0) for a given /i, the frequencies, damping 
rates and frequency ratio are greater than those of a major one. However, the results obtained 
for the aforementioned quantities in the major elliptical emerged loop are slightly greater than 
those of a circular-arc semi-emerged loop (e = A = 0). 

v) The ratio of the oscillation frequency to the damping rate, w/|7|, in minor /major elliptical 
and circular-arc semi-emerged loops, is not affected by making changes in density stratification 
parameter, ellipticity and stage of emergence of the loop. 
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Figure 1: Frequency of the fundamental kink {m = 1) mode and its damping rate as well as the 
ratio of the oscillation frequency to the damping rate as a function of the stratification parameter 
/X for a circle-arc flux tube (e = 0) at various stages of emergence A = (sohd), 0.75 (dashed) 
and -0.75 (dash-dotted). The loop parameters arc L = 10^ km, R/L = 0.01, l/R = 0.02, 
Pcx/Pin = 0.1, Pin = 2 X lO^^^gr cm^'^ and B = 100 G. Both frequencies and damping rates are 
in units of the interior Alfven frequency, loa^^ = 0.02 rad s~^. 
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Figure 3: Ratio of the frequencies ui2/^i of the first-overtone and its fundamental kink {m = 1) 
mode versus /v, for a circle-arc flux tube (e = 0) at various stages of emergence A = (solid), 
0.75 (dashed) and —0.75 (dash-dotted). Auxiliary parameters as in Fig. 1. 
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Figure 4: Frequency of the fundamental kink (m = 1) mode and its damping rate as well as 
the ratio of the oscillation frequency to the damping rate as a function of the stratification 
parameter fj, for a minor scmi-cmcrgcd loop (A = 0) with different ellipticity parameters e = 0.0 
(solid), 0.4 (dashed) and 0.6 (dash-dotted). Auxiliary parameters as in Fig. 1. 
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Figure 6: Ratio of the frequencies io^j^x of the first-overtone and its fundamental kink (m = 1) 
mode versus /v, for a minor scmi-cmcrgcd loop (A = 0) with different ellipticity parameters e = 0.0 
(solid), 0.4 (dashed) and 0.6 (dash-dotted). Auxiliary parameters as in Fig. 1. 
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minor case, e =0.6, A.=0.75 
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Figure 7: Frequency of the fundamental kink (m = 1) mode and its damping rate as well as the 
ratio of the oscillation frequency to the damping rate as a function of the stratification parameter 
12 for circular-arc semi-emerged (solid), minor (dash-dotted) and major (dashed) elliptical loops. 
Auxiliary parameters as in Fig. 1. 
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major case, e =0.6, X=0.75 
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Figure 9: Ratio of the frequencies uj2/iOi of the first-overtone and its fundamental kink (m = 1) 
mode versus /i for circular-arc semi-emerged (solid), minor (dash-dotted) and major (dashed) 
elliptical loops. Auxiliary parameters as in Fig. 1. 
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